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Abstract: Detached-eddy simulation (DES) based on the v 2 -f Reynolds-averaged Navier-Stokes 
(RANS) model is developed and tested. The v 2 -f model incorporates the anisotropy of near- wall 
turbulence, which is absent in other RANS models commonly used in the DES community. The 
v 2 -f RANS model is modified in order that the proposed u 2 -/-based DES formulation reduces 
to a transport equation for the subgrid-scale kinetic energy in isotropic turbulence. First, three 
coefficients in the elliptic relaxation equation are modified, which is tested in channel flows with 
friction Reynolds number up to 2000. Then, the proposed v 2 -f DES model formulation is derived. 
The constant, Cdes? required in the DES formulation was calibrated by simulating both decaying 
and statistically-steady isotropic turbulence. After Cdes was calibrated, the v 2 -f DES formulation 
is tested for flow around a circular cylinder at a Reynolds number of 3900, in which case turbulence 
develops after separation. Simulations indicate that this model represents the turbulent wake nearly 
as accurately as the dynamic Smagorinsky model. Spalart-Allmaras-based DES is also included in 
the cylinder flow simulation for comparison. 
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1 Introduction 

Detached-eddy simulation (DES) [1, 2, 3] is a hybrid RANS/LES approach that performs RANS in attached 
regions and LES in detached regions using a single model. In the LES region, the length scale of the model 
is set proportional to the grid size A. The RANS model thereby becomes an LES model. 

Although the Spalart-Allmaras (SA) model has been widely used for DES [3], its near- wall damping does 
not distinguish between velocity components. In contrast, the v 2 -f formulation [4] models the suppression of 
wall normal velocity fluctuation caused by non-local pressure-strain effects. This anisotropy has been shown 
to improve prediction of separation and reattachment [5, 6]. Among the few works where both v 2 -f and SA 
RANS were compared, Iaccarino et al. [6] predicted the pressure coefficient C p with v 2 -f more accurately 
than with SA for flow over a 2D hump with steady-suction flow control. Const antinescu et al. [7] achieved an 
accurate friction coefficient Cf on a sphere after laminar separation with v 2 -f but an inaccurate prediction 
with the k-w, fc-e, and SA RANS models. 

The u 2 -/-based DES model developed here incorporates more flow physics than SA-based DES. The SA 
model uses the minimum distance to the wall as the turbulence length scale which is not necessarily accurate 
at or near separation. The v 2 -f model, on the other hand, computes a length scale based on flow properties 
(namely, the kinetic energy fc, the dissipation rate e, and the kinematic viscosity v). In addition, the present 
v 2 -f DES model gives a transport equation for subgrid-scale (sgs) kinetic energy which is less empirical than 
the sgs viscosity transport equation used in SA-based DES. 

A brief summary of the proposed v 2 -f DES methodology is as follows. The v 2 -f RANS model [4] has 
transport equations for fc, e, and v 2 together with an elliptic relaxation equation for a function /. These 
equations contain length and time scales L and T. The choice between RANS and LES modes is made by 
setting these scales appropriately. When the grid is fine enough to capture large-scale turbulent eddies, in 
particular, when Cdes A < k 3 / 2 /e, LES mode is selected. Otherwise, RANS mode is selected. As suggested 
by Spalart et al. [1], the coefficient Cdes is chosen to match the correct energy spectrum in isotropic 
turbulence. In LES mode, the time scale is set to T = Cdes A/ where k now represents the sgs kinetic 
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energy. Although the k equation is sufficient for purely LES simulations, the other three equations are 
required for the RANS mode. The v 2 -f model [4] is modified so that the entire set of equations reduces in 
LES mode to a transport equation for the sgs k. In particular, coefficients in the elliptic relaxation equation 
are modified so that v 2 becomes statistically (2/3 ) k in the limit of LES for isotropic turbulence. 

The proposed v 2 -f DES model is described in section 2. In section 3, results for this model are compared 
with the dynamic Smagorinsky model [8] and SA-based delayed DES [2] using the same flow solver and grid 
for flow past a cylinder at Reynolds number based on the diameter Rep = 3900. This cylinder flow is a good 
test case with grids appropriate for non-dissipative schemes. This test is not chosen to determine which DES 
model is better than others, since it is expected that all DES models will behave like the Smagorinksy model 
for this flow. Conclusions are presented in section 4. 


2 v 2 -f DES Model 


The DES model proposed is based on Lien et al.’s version of the v 2 -f RANS model [4]. The quantity 
v 2 should be interpreted as a velocity scalar, not a Reynolds stress tensor component, and the function / 
makes v 2 behave like the wall normal component of the Reynolds stress (u' n u' n ) [5]. Another meaning that we 
believe should be imposed on this scalar is that in isotropic turbulence v 2 should be the average of the normal 
components of the Reynolds stress. This assumption allows us to formulate a DES model that reduces to a 
transport equation for sgs -k in isotropic turbulence. This derivation starts from the v 2 -f RANS Eqs.(l)-(5). 
To reach the final formulation (Eqs.(18)-(22)), the formulations of sgs -k LES (Eqs.(12)-(14)) and k-e DES 
(Eqs.(15)-(17)) are derived on the way. 

The v 2 -f RANS model of Lien et al. [4] is 
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Note that the realizability constraints of Durbin [9] are imposed to both scales [10]. The turbulence produc- 
tion V = 2v t \Sij\ 2 , and the strain rate tensor Sij = \ ( djUi + diUj). The coefficients are given by 


c„ = 0.22, Cel = 1.4 (l + 0.045 ^J) , c e2 = 1.9, a e = 1.3, 

c\ = 1.4, C 2 = 0.3, ct = 6, cl = 0.23, c ^ = 70. (7) 

Wall boundary conditions are 

kw = 0, € w = vd n k\ w , v 2 ^ : 0, f w - 0, (8) 

where d n is the wall-normal gradient. Note that in Eq.(5), we have chosen to keep cl outside of the definition 

of Trans for a later purpose, unlike Lien et al. [4]. 
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To begin with, we modify the three coefficients in the elliptic relaxation equation (5) as follows, 


C2 = C2,0 + |(2/3 — C2,o) 


- tanh jc 2 ,i(u 2 /& - c 2 ,2)| , c 2 ,o = 0.3, c 2 ,i = 50, c 2 ,2 = 0.55 
cl = 0.115, c jj = 140 


(9) 

(10) 


The value of C 2 ,o is the unmodified value of C 2 . The other parameters C 2 ,i and 02,2 are chosen for the 
modified C 2 — > 2/3 when v 2 /k > 0.6, as shown in Fig. 1. This increases the rapid distribution term (the 
last term in Eq. (5)) so that it reduces to the isotropization of production (IP) model [11] when C 2 = 2/3. 
The parameters C 2 ,i and C 2,2 are not optimized. Effects of these parameters on a turbulence wake behind a 
cylinder are discussed in section 3. The modification of C 2 allows v 2 to be statistically 2/3 k when this model 
is used as an sgs model for isotropic turbulence, i.e., (v 2 ) = (2/3) (k). 



Figure 1: The modified C 2 of the present v 2 -f RANS model 

To explain how the v 2 equation (4) becomes 2/3 of the k equation (2), let us approximate / as / where 
/ is defined by 
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This approximation is valid for y + > 10 in channel flow, as shown in Fig. 2. The channel flow simulation is 
discussed in detail later. Substituting this approximation into Eq. (4), the right-hand side of the v 2 equation 
(4) becomes 2/3 of the k equation (2) because c 2 — > 2/3. It was found that the modified c 2 (Eq. (9)) 
required cl and c v to also be modified to yield acceptable results for channel flow. The value of the product 
clc^ — 16.1 is the same in order to keep the Laplacian term near the wall the same as in the unmodified 
form and to reduce it away from wall. 

As briefly mentioned above, to validate the modified coefficients of the elliptic relaxation equation, fully 
developed turbulent channel flows with friction Reynolds numbers Re r up to 2000 are simulated in RANS 
mode. Figure 3 shows that this modification maintains the performance of the unmodified RANS model 
with a slightly better prediction of v 2 near the center. 

In the DES approach, a RANS model switches to an LES model based on a comparison of length scales. 
The LES length scale is determined by the grid size A, so only a time scale equation is needed for an sgs model. 
We would like the v 2 -f DES model to reduce to an sgs -k equation model similar to Yoshizawa’s model [12] 

given that for isotropic turbulence (v 2 ) becomes the average of the sgs normal stress i.e., (v 2 / = (2/3) (k). 
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Figure 2: The function / + ( — ) and its approximation / + = 5 v 2+ e + /k +2 J cC 2 V + /k Jr i.e., Eq.(ll) (- -) in 
the present v 2 -f RANS for the channel. The profiles of Re r = 550, Re r = 950 and Re r = 2000 are shifted 
upward along the vertical axis by 0.01, 0.02 and 0.03, respectively. The plus superscript indicates wall units. 






Figure 3: Comparison of DNS’s (dots) and RANS with the unmodified [4] ( — ) and the present v 2 -f model 
(— • — ) in turbulent channel flows. The plus superscript indicates wall units. 
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An sgs -k LES model is 
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When the flow is in equilibrium i.e., V = zees, a Smagorinsky-like eddy viscosity v t oc \Sij\A 2 is obtained. 

As an intermediate step towards the v 2 -f DES model, a simplified k-e based DES model (without damping 
functions) is presented. Incorporating the dissipation equation into the DES context yields 
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If the RANS length scale is selected, then Tees = k/e and cdes = e as expected. Otherwise, Tees = &Ales 
which is the same as Tles in Eq.(12), and Eqs.(15)-(17) reduce to Eqs.(12)-(14). 

Finally, to present the full v 2 -f DES model, the v 2 and / equations are included in the complete set of 
equations: 
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if k 3 / 2 /e < Lees, then Lees = Leans, Tees = ^ranSj anc i 6 des = e (23) 

otherwise, Lees = CdesA, Tees = CdesA/v^/c, and cdes = ^ 3 ^ 2 /(C f DEsA), (24) 


where Trans an d Trans ar *e defined in Eq.(6). Note that (2/3 ) k appears in Eqs.(12) and (15) but v 2 appears 
in Eq.(18). For an LES simulation of isotropic turbulence, u 2 , if interpreted as the sgs average normal stress, 
should be statistically (2/3) k. The modified c 2 in Eq.(9) provides this behavior. Simulations indicate that 

the unmodified c 2 = 0.3 yields / (k) — 0.4 for isotropic turbulence. 

Note that the Kolmogorov length scale present in Eq. (6) is not considered in the length comparison 
in Eqs. (23) and (24) because it is expected that, near the RANS and LES transition (or away from the 
wall), the RANS length scale will be k 3 ! 2 /e. However, discontinuity in length or time scale may occur if the 
Kolmogorov length scale happens to be the RANS scale at the RANS/LES transition. This would occur 
in low-Te flows. Another way to avoid such a discontinuity is maybe to use Trans for the length scale 
comparison. This alternative is not considered here because it turns off the realizability constraints in Eq. 
(6) causing the stagnation point anomaly [9] near the front stagnation point on a cylinder. Flow simulation 
over a cylinder is discussed in section 3. 

The proposed DES model is implemented in the Stanford incompressible Navier-Stokes solver CDP 
(v2.3). This code is based on finite- volume spatial discretization with second-order accuracy and a second- 
order implicit fractional step method [13, 14]. CDP employs a novel collocated formulation to conserve mass, 
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momentum and approximately kinetic energy (in the inviscid limit) on a hybrid unstructured mesh. The 
second-order central difference scheme is used for both the convection and diffusion terms of the momentum 
equation. The scalar transport equations (19)- (21) are time advanced implicitly. To achieve numerical 
stabilization for the scalar equations, which are significantly controlled by the source terms, the convective 
terms in these equations are discretized with the lst-order upwind scheme. The diffusion term in each scalar 
equation is discretized with the central scheme. In the test of isotropic turbulence, when this upwind scheme 
is used for the scalar convection term only, it yields the identical energy spectrum compared to use of the 
central scheme for all the terms. An upwind-biased scheme is never used for the momentum equation in this 
study. 

To calibrate Cdes? decaying isotropic turbulence is considered. The initial field uses 512 3 DNS data of 
Wray [15] at Re \ = 105 sampled down to 32 s . The computational cube has sides Lbo X = 27r with periodic 
boundary conditions in the three directions. Quantities in the DES model are initialized by running for some 
time with a frozen initial velocity field. The flow is allowed to run about 3.7 and 10.7 large-eddy turnover 
times at which point Re \ = 65 and Re \ = 61, respectively. The large-eddy turnover timescale is determined 
by the longitudinal length scale L\\ and the r.m.s. velocity fluctuation u rms computed from the DNS data 
at Re a = 105 (the initial state). 

Figs. 4 (a) and (b) show the energy spectrum with three values of Cdes- The v 2 -f DES simulation with 
Cdes = 0-8 agrees very well with DNS [15] and experiments [16] up to the cut-off wave number at the later 
instant when Re \ = 61 (Fig. 4 (b)). v 2 -f DES gives almost exactly the same energy spectrum as produced 
by the dynamic Smagorinsky model (Figs. 4 (c) and (d)). In Figs. 4 (e) and (f), all three formulations, sgs -k 
LES, k-e DES, and v 2 -f DES, give the same energy spectrum with the same Cdes- This indicates that the 
v 2 -f DES formulation indeed reduces to the others. The coefficient Cdes = 0-8 is used for the rest of the 
paper. 

Because it allows a higher Reynolds number and a wider inertial range, forced isotropic turbulence at 
Re \ = 98 is also considered using the stochastic forcing of Eswaran and Pope [17]. The computational cube 
has the same size Lbo X = 27r as the decaying case with the same periodic conditions. The radius of the 
sphere of forced wavenumbers is Kp = \/8, giving a total 92 forced modes. The forcing amplitude and the 
forcing autocorrelation time scale are chosen as a = 0.3572 and T L = 0.4312, respectively, following the 
nomenclature of Ref. [17]. Initial conditions are not relevant, since the flow is driven to a statistical steady 
state. The grid is 128 3 for the present DNS, and two grids of 16 3 and 32 s are used for both LES with the 
dynamic Smagorinsky model and v 2 -f DES. Statistics were obtained by averaging over samples collected 
over 30 eddy turnover times in the steady state which is reached after 10 eddy turnover times from the start. 
The velocity-derivative skewness 5 = 1/3 (( diU { ) 3 ) / (( djUj ) 2 ) 3 / 2 has values (^ 0.45) typical of isotropic 

turbulent at Re \ ~ 100 [18]. The v 2 -f DES simulation gives almost exactly the same energy spectrum as 
produced by the dynamic Smagorinsky model (Fig. 5). Forced isotropic turbulence had not been used for 
the Cdes calibration in the DES community. This simulation shows that the forced case is an alternative 
for this calibration. 
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(a) v 2 -f DES with three different Cdes at Re\ = 65 



(b) v 2 -f DES with three different Cdes at Re\ = 61 



(c) LES and v 2 -f DES with C D es = 0.8 at Re x = 65 



(d) LES and v 2 -f DES with Cdes — 0.8 at Re\ = 61 




(e) sgs -k LES, k-e DES and v 2 -f DES with Cdes — 0.8 (f) sgs -k LES, k-e DES and v 2 -f DES with Cdes — 0.8 


Figure 4: Energy spectrum of decaying isotropic turbulence. The proposed v 2 -f DES model is simulated 
with three different values of Cdes on 32 3 in (a) and (b). The curve labeled “-5/3” shows E(k) = Ce 2 / 3 ^ -5 / 3 
with the Kolmogorov constant C = 1.5. “DynSmag” indicates LES with the dynamic Smagorinsky model. 
The coefficient Cdes = 0-8 is used for sgs -k LES, k-e DES and v 2 -f DES on 32 s in (e) and (f). 
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Figure 5: Energy spectrum of forced isotropic turbulence of Re \ = 98. Shown are the present simulations 
of DNS, LES with the dynamic Smagorinsky model, and v 2 -f DES with Goes = 0.8 on the 16 3 grid (left) 
and the 32 s grid (right). 


3 Results: Flow over a circular cylinder 

To test the LES mode of the proposed v 2 -f DES model, flow over a circular cylinder is chosen. The Reynolds 
number based on diameter D is Rep = 3900. The DES model coefficient Cdes 5 the only coefficient which 
has a degree of freedom, is determined in isotropic turbulence. Therefore, it is imperative to test the DES 
model in anisotropic turbulence. This flow over a circular cylinder is a good test case, because (1) turbulence 
in the wake is anisotropic, and (2) this flow has been extensively studied both experimentally [19, 20] and 
numerically [14, 21, 19, 20]. Kravchenko and Moin [21] used a Galerkin B-spline method with the dynamic 
Smagorinsky model to simulate this flow with good agreement of mean flow quantities with experimental data 
available at the time. Mahesh et al. [14] simulated this flow with the dynamic Smagorinsky model to validate 
the numerical method used in the CDP code. Dong et al. [19] performed DNS using a spectral element 
scheme on an unstructured grid and compared it with their PIV experiments. Recently, Parnaudeau et al. 
[20] performed an extensive study in order to address the lack of consensus in the literature for turbulence 
statistics immediately behind the cylinder. Parnaudeau et al. [20] simulated this flow with the dynamic 
Smagorinsky model using an immersed boundary method on uniform grids, and compared the computations 
with their PIV and hot-wire experiments. 

The cylinder flow at Rep = 3900 allows us to test DES as an LES model with grids appropriate for 
non-dissipative schemes. Strelets [22] and Travin et al. [23] performed SA-based DES for a cylinder at 
Rcb = 50000. We attempted this higher - Re case first but discovered that the resolution used in the previous 
studies with dissipative upwind-biased schemes was inadequate for the present second-order non-dissipative 
code. Strelets used a hybrid central/upwind approximation [22], and Travin et al. used a fifth-order upwind 
scheme [23]. The resolution in Strelets [22] and Travin et al. [23] is even coarser than that in previous LES 
computations at Rev = 3900 [24, 21, 14, 20]; Refs. [22, 23] used about 0.7 million points, whereas 1.3 - 
4.4 million points are used in Refs. [24, 21, 14, 20]. Mittal and Moin [24] concluded that even a 5th-order 
upwind scheme can degrade LES computations at Rep = 3900. Simulation of the Rep = 50000 flow is still 
challenging with non-dissipative schemes. Because both Re& = 3900 and 50000 are subcritical (turbulence 
develops after separation), the Rep = 3900 flow is a good test case with grids appropriate for non-dissipative 
schemes. This test is not chosen to determine which DES model is better than others, since it is expected 
that all DES models will behave like the Smagorinksy model for this flow. 

The present computation was performed on a domain whose inflow and outflow surfaces were 30 D up- 
stream and 35 D downstream from the center of the cylinder, following the grid of Mahesh et al. [14] . The 
domain height was 50 D and spanwise extent was ttD. The size of the first grid adjacent to the cylinder 
is 0.0025D radially and 0.01D (or 0.57°) in the azimuthal direction 0. The quadrilateral elements are ap- 
proximately 0.04D x 0.04D at a distance of 2D from the cylinder center in the wake. The current grid 
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has slightly higher resolution on the wall than grids in previous LES studies [21, 14, 20] and less resolution 
than the DNS study [19], particularly in the azimuthal 0 and spanwise z directions. There are twelve grid 
points inside the boundary layer at separation. Fig. 6 shows the domain and the grid. The flow field is 
initially uniform. Without any imposed perturbation, simulation generates 3D fluctuations due to biased 
round-off error. Note that this flow is naturally unstable to perturbation (see the Hopf bifurcation theory 
[25]). A statistically steady state is obtained after about 200 time units. Statistics were accumulated over 
approximately 30 vortex shedding cycles (140 time units) and over the spanwise direction. The time step is 
A t = O.OOID/Uoo. Table 1 summarizes the computation parameters in the previous and current studies. 



Figure 6: Computation domain (left) and grid around the cylinder (right). 


Study 

Kravchenko et al. [21] 

Mahesh et al. [14] 

Parnaudeau et al. [20] 

Dong et al. [19] 

Present 

Domain [D] 

(R d ,L z ) = (30,7 r) 

Lx,y,z — (75, 50, 7f) 

L x ,y,z = (20, 20, tt) 

L x ,y,z = (40, 18, 7r) 

Lx,y,z = (75, 50, 7r) 

Grid 

n r ,e,z = (205, 185,48) 

n z - 32* 

n x , y ’ z = (961, 960, 48)* 

n z - 128* 

n z = 64* 

A r, el wall [10“ 2 D] 

(0.25, 1.7) 

(0.25, 1.2) § 

(2.1, 2.1) 

(0.36, 0.77)^ 

(0.25, 1.0) 11 

A t [10 -Zd/Uoo] 

2.5 

1 

3 

2 

1 

Simulation type 

LES 

LES 

LES 

DNS 

LES, DES 


Table 1: Computation parameters used for the current study and previous studies. Note *: unstructured 

quadrilateral grid with resolution of 0.04D for a distance of 1.5D in the wake, uniform mesh, unstructured 
triangular grids, minimum values on the cylinder, calculated from their values in viscous wall units, H; 
uniform on the cylinder. 

Three turbulence models are used for the present computations: dynamic Smagorinsky LES [8], SA-based 
delayed DES (SA-DDES) [2] and the proposed v 2 -f DES model. The SA equation without the trip term 
[26], called SA-noft2 in the turbulence community, is selected here because this version can use lower values 
of v at freestream boundary [27]. The low- Reynolds number correction of Spalart et al. [2] is required for 
SA-DDES to compensate for unexpected activation of the wall functions possibly due to either low Reynolds 
number or very fine grids. This SA-DDES was also implemented in CDP. Using Cdes = 0-65 for SA-DDES 
gives the correct energy spectrum in isotropic turbulence (data not shown). A = max(A^) in all cases. Initial 
and boundary conditions are listed in Table 2. 

Instantaneous velocity and vorticity fields in the wake of a circular cylinder at Ren = 3900 are shown in 
Figs. 7 and 8, respectively. Only the results of the v 2 -f DES simulation are shown in these instantaneous 
figures, and statistics from three models are compared later. Fig. 7 (a) shows unsteady recirculation region 
behind the cylinder. Alternating regions of positive and negative cross-stream velocity are shown in Fig. 
7 (b). This is related to Karman vortices shown in Fig. 8. Fig. 7 (c) clearly shows fluctuating spanwise 
velocity which indicates three dimensional flow structures in the wake. Figs. 7 and 8 are conformable to the 
LES results of Kravchenko et al., i.e., Figs. 2-5 of Ref. [21], which show the similar large structures. 

Mean flow fields are obtained from the simulation over approximately 30 vortex shedding cycles and over 
spanwise direction in the near wake of the cylinder. Fig. 9 shows the results of v 2 -f DES. The velocity deficit 
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Variables 

All 

v z -f DES 

SA-DDES 

Ui 

k 

v 2 

e 

f 

V 

Initial Conditions 

uniform Uoo 

10 -a ' 

(2/3) x 10- 

^ 10 

none 

hv 

Wall Boundary 

0 

0 

0 

vd z k 

0 

0 

Inlet Boundary 

Uoo 

0 

0 

zero gradient 

zero gradient 

0.5 v 

Outlet Boundary 



convective outlet 


zero gradient 

convective outlet 

Slip Boundary 

zero gradient 

Z planes 

periodic 


Table 2: Initial and Boundary conditions in the present computations. The gradient d n is in the wall normal 
direction. 



Figure 7: Instantaneous velocity fields on the spanwise mid-plane (y = 0) in the wake of a circular cylinder 
at Rev = 3900 from the v 2 -f DES simulation. Shown are the normalized streamwise velocity uju^ (a), the 
normalized cross-flow velocity v/u^ (b), and the normalized spanwise velocity w/u^ (c). The thick black 
lines indicate the zero value. The number of color contour levels is 52 between —1.5 (blue) and 1.5 (red). 
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Figure 8: Instantaneous vorticity magnitude loD/uqq on a x-y plane with 16 contours from loD/u^ = 0.5 to 
c jD/uoq = 10.0 in the simulation of v 2 -f DES. 


in the wake is well shown in Fig. 9 (a) with the negative bubble near the cylinder. Separating shear layer is 
clearly shown in the spanwise vorticity field, Fig. 9 (b). The Reynolds shear stress (u'v') is anti-symmetric 
with respect to the centerline (see Fig. 9 (c)). (u'v') has has two dominant peaks of | (u'v') | = 0.119 at 
x/D = 2.08 and two secondary peaks of | (u'v') | = 0.040 at x/D = 1.37. The r.m.s. streamwise velocity 
u rms in Fig. 9 (d) shows strong gradient in the separating shear layers and two maxima associated with 
the vortex formation in Fig. 9 (b). These mean flow structures are qualitatively similar to Parnaudeau et 
al.[20]’s experimental and LES data, and Dong et al. [19] ’s experimental and DNS data. Fig. 10 accumulates 
present simulation results with fewer contour fines for direct comparison among the three models. In general, 
all three models provide the similar mean flow fields with subtle differences which are emphasized by the 
same level of the contour fines. The simulation without any turbulence model is not included in Fig. 10 
because its results are simply not correct (see Table 3). 

Integral flow quantities are well predicted in the present computations with three models, as fisted in 
Table 3. The drag coefficient Cd = 1-00 of v 2 -f DES matches well with previous experiment data [21] and 
computations [14]. This coefficient is overestimated without any turbulence model. The value Cd = 0.965 
predicted by SA-DDES is acceptable. The base pressure coefficient C p ^ = —0.928 of v 2 -f DES is similar to 
that of Kravchenko et al. [21]. SA slightly underestimates C p ^> and, without any turbulence model, C p ^ is 
simply wrong. The separation point 0 sep = 87.1 of the present dynamic-Smagorinsky LES and v 2 -f DES is 
in the range of the experimental data set and similar to previous LES results [21, 14, 20]. Previous studies 
show a broad range of the recirculation length Lr ; this value is sensitive to experimental conditions [20], and 
previous LES slightly underestimated this value. The current value Lr = 1.44 of v 2 -f DES is in the range 
of the experimental data. Without a turbulence model, the recirculation length is completely wrong. The 
Strouhal number St is the least sensitive quantity in Table. 3. The current LES, v 2 -f DES and even the 
simulation with no model predict the shedding frequency very accurately. SA-DDES slightly overestimates 
St. The minimum streamwise velocity U m i n predicted by all the present computations is within the range 
of the previous data. 

Fig. 11 shows mean velocity in the wake. v 2 -f DES predicts the mean streamwise velocity on the 
centerline observed in previous studies [21, 20] with slightly better prediction downstream 5 < x/D < 10 
than the other present computations (Fig. 11(a)). All current LES and DES’s reproduce the mean streamwise 
and cross-flow velocity at three locations in the wake as shown in Fig. ll(b,c). For the mean cross-flow 
velocity V at x/D = 1.54, the present LES and v 2 -f DES match closer with the experiment of Parnaudeau 
et al. [20] than does their LES. 

Fig. 12 shows Reynolds stress profiles in the near wake. Overall, all three models predict the three 
components quite well, compared with LES and experimental data of Parnaudeau et al. [20]. At x/D = 1.54, 
v 2 -f DES and the present dynamic-Smagorinsky LES are slightly better than SA-DDES, particularly for 
( u'u ') and (u'v'). For ( u'u '), v 2 -f DES and the present dynamic-Smagorinsky LES match closer with the 
experiment of Parnaudeau et al. [20] than does their LES, similar to V in Fig. 11 (c). Such better prediction 
in the present simulation is probably due to better resolution very near the cylinder; the present grid is finer 
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Figure 9: Mean flow fields in the simulation of v 2 -f DES: (a) the normalized streamwise velocity with 
contour of f//^oo|min = —0.3 and lAU/u^l = 0.1; (b) the normalized spanwise vorticity with contour 
of | (u z ) EV^oolmin = 1 and | A (lj z ) D/u ^ | = 1; (c) the normalized Reynolds stress with contour of 
| (u'v') = 0.018 and |A (u'v') /u^l = 0.009; and (d) the normalized r.m.s. streamwise velocity 

with contour of ^rms/^oolmin = 0.04 and lA'Urms/'Uool = 0.04. In all plots, dashed lines indicate negative 
contour. 



Study 

Cd 

Cp,b 

$sep 

Lr/d 

St 

Uni in / Uoo 


Data set* 

0.99T0.05 

-0.88T0.05 

86.0±2 

1.4T0.1 

0.215T0.005 

-0.24i0.01 

Exp. 

Parnaudeau et al. [20] 



88 

1.51 

0.208T0.002 

-0.34 


Dong et al. [19] 




1.47 


-0.252 

DNS 

Dong et al. [19] 


-0.88 


1.36 

0.20 

-0.291 


Kravchenko et al. [21] 

1.04 

-0.94 

88.0 

1.35 

0.210 

-0.37 

LES 

Mahesh et al. [14] 

1.00 


87.6 

1.35 

0.218 

-0.31 

Parnaudeau et al. [20] 



87.6 

1.35 

0.208 

-0.31 


Present 

0.994 

-0.881 

87.1 

1.48 

0.214 

-0.320 

DES 

Present, v z -f DES 

1.00 

-0.928 

87.1 

1.44 

0.214 

-0.304 

Present, SA-DDES 

0.965 

-0.969 

88.3 

1.37 

0.221 

-0.283 


Present, w/o turb. model 

1.18 

-1.16 

89.4 

0.841 

0.214 

-0.264 


Table 3: Mean flow quantities of the flow around the circular cylinder at Ren = 3900: the drag coefficient 
Cd , the base pressure coefficient C Pi b, the separation angle from the front # sep , the recirculation length L#, 
the Strouhal number St , and the minimum streamwise velocity U m j n on the centerline. Note *: experimental 
data in Table II of Kravchenko et al. [21] 
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Figure 10: Mean flow fields in the simulation of v 2 -f DES (-), dynamic-Smagorinsky LES (-), and SA- 
DDES (-): (a) the normalized streamwise velocity with contour of U/u 0 o| m in = — 1 0.2 and {AU/u^l = 0.2; 
(b) the normalized spanwise vorticity with contour of | (u z ) -D/^oolmin = 2 and |A (uj z ) D/u^l = 4; (c) the 
normalized Reynolds stress with contour of | ( u'v ') /u^ | m in = 0.018 and |A ( u'v' ) /u^ \ = 0.036; and (d) the 
normalized r.m.s. streamwise velocity with contour of 'Urms/'^oolmm = 0-08 and {Au^s/u^l = 0.12. In all 
plots, dashed lines indicate negative contour. 
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(a) U/Uoo (b) U/zioo (c) V/ Uoo 

Figure 11: Mean velocity profiles in the wake of the circular cylinder at Rer> = 3900. Shown are (a) 
normalized streamwise velocity U/u^ on the centerline, and (b) U/u oo and (c) normalized cross-flow velocity 
V/u^ at three locations in the wake of a circular cylinder at Ren = 3900. In (b) and (c), profiles at 
x/D = 1.54 and x/D = 2.02 are vertically shifted down by 1 and 2, respectively. 



(a) (u'u f ) /ulv (b) ( u'v ') /ulc (c) (v'v') / u ^ 

Figure 12: Normalized Reynolds stress profiles at three locations in the near wake of a circular cylinder at 
Reo = 3900. In (a), profiles at x/D = 1.54 and x/D = 2.02 are vertically shifted down by 0.3 and 0.6, 

respectively. In (b), profiles at x/D = 1.54 and x/D = 2.02 are vertically shifted down by 0.2 and 0.4, 

respectively. In (c), profiles at x/D = 1.54 and x/D = 2.02 are vertically shifted down by 0.4 and 0.8, 

respectively. See Fig. 11 for labels. 
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than that of Parnaudeau et al. [20] from the cylinder to x/D ~ 1.2 in the wake. 

Compared to v 2 -f DES and dynamic-Smagorinsky LES, SA-DDES predicts rather smooth profiles of V, 
{ u'u '), and ( u'v ') in the near wake. This is related to large mean eddy viscosity (v t ) produced by SA-DDES, 
as shown in Fig. 13. v 2 -f DES and dynamic-Smagorinsky LES generate (u t ) ~ v in the near wake around 
1 < x/D < 2, whereas (v t } ~ from SA-DDES. Interestingly, both v 2 -f DES and SA-DDES have (v t ) 
peaks around x/D ~ 1.5 close to the centerline in the near wake, whereas the dynamic Smagorinsky has two 
symmetric peaks at x/D ~ 2.5 and y/D ~ ±0.5. Although v 2 -f DES decreases (v t ) downstream x/D > 2, 
small (v t ) of v 2 -f DES does not significantly affect mean flow field in the far wake (see Fig. 15). 





(c) SA-DDES 

Figure 13: Contours of the normalized mean eddy viscosity (u t ) /v. In all plots, (v t ) /^| m in = 0-2 and 
| A (v t ) / v\ = 0.2. 

Fig. 14 shows mean model quantities of the v 2 -f DES model in the flow over the cylinder. Overall, a 
similar structure is observed in (fc), and (e); they are symmetric with respect to the centerline with 

peaks around x/D ~ 1.5 near the centerline, (cdes) l ar g er than (e), which indicates that the LES mode is 
active in the wake and that the LES dissipation is larger than the RANS dissipation, as expected. The ratio 
/ (k) increases from almost zero near the wall to 0.42 around x/D ~ 1 and slowly decreases as the eddy 
viscosity decreases. Since v 2 /k does not reach 2/3 in this simulation, C 2 — C 2 ,o in the whole domain. This 
phenomena is rather unexpected, although ^u 2 ^ / (k) = 2/3 is the asymptotic value in the limit of isotropic 

LES. In order to make ^ 2 ^ / (&) — > 2/3 in the wake, the parameters c 2 ,i and c 2 ,2 in Eq. (9) probably need to 

be modified, so C 2 changes more smoothly between C 2 ,o to 2/3. As mentioned in section 2, these parameters 
are not optimized. In Fig. 14 (e), max((fc) , le — 6) is used for (k) to avoid 0/0. The mean function (/) 

has peaks similar to those of ^ 2 ^> which visually shows that / suppresses v 2 near the wall and produces v 2 
away from the wall. 

Mean flow quantities in the far wake are shown in Fig. 15. Only the data of Kravchenko et al.[21] is 
included as the reference, because Parnaudeau et al.[20] did not include far wake data. Data of Kravchenko 
et al. [21] has good agreement with the experiment of Ong and Wallace [28]. Similar to the near wake, all 
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(e) (v 2 ^ / (k) 


(f) </> D/u a 


Figure 14: Normalized mean model quantities of the v 2 -f DES model, (a) |A (fc) = 0.001. (b) 

|A (f) A4>l = 0-005. (c) |A (e DES ) D/ulcl = 0.003. (d) |A (e) D/v.%, \ = 0.003. (c) |A (v 2 ) /(k)\= 0.05. 
(f) |A (/) U/iXool = 0.6. 
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three turbulence models give the acceptable profiles. v 2 -f DES has slightly better agreement in (u'v') than 
the others, whereas dynamic Smagorinsky is slightly better for ( u'u '). Based on both ( u'u' ) and (v'v'), this 
far-wake flow is still far from isotropic turbulence. 




y/D y/D 


(c) (u'u') /ule (d) (v'v') /ulc 

Figure 15: Mean streamwise velocity (a), Reynolds shear stress (u'v') (b), and normal stresses (u'u') (c) 
and (v'v') (d), at three locations in the far wake of a circular cylinder at Reo = 3900. In (a), profiles at 

x/D = 7.00 and x/D = 10.0 are vertically shifted down by 0.3 and 0.6, respectively. In (b), profiles at 

x/D = 7.00 and x/D = 10.0 are vertically shifted down by 0.01 and 0.02, respectively. In (c), profiles at 

x/D = 7.00 and x/D — 10.0 are vertically shifted down by 0.05 and 0.1, respectively. In (d), profiles at 

x/D = 7.00 and x/D = 10.0 are vertically shifted down by 0.1 and 0.2, respectively. See Fig. 11 for labels. 
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4 Conclusions and Future Work 


DES based on the v 2 -f model was proposed, implemented in the unstructured incompressible code CDP, 
and tested for isotropic turbulence and flow around a circular cylinder. Three coefficients in the elliptic 
relaxation equation of the original RANS model are modified in order that v 2 is statistically (2/3) k in the 
limit of isotropic turbulence. This allows the v 2 -f DES formulation to reduce to sgs -k LES in this limit. The 
DES coefficient Cdes determined from isotropic turbulence is Cdes = 0-8- 

Flow around a circular cylinder at Rep = 3900 is simulated with the proposed v 2 -f DES along with SA- 
based DDES and the dynamic Smagorinsky model. Since at this Reynolds number turbulence occurs after 
separation, this case tests only the LES mode. v 2 -f DES reproduces not only instantaneous large structures 
in the wake but also mean flow fields as observed in previous experimental and numerical studies. Dynamic- 
Smagorinsky LES and SA-DDES are performed for direct comparison among the three models. Overall, 
all three turbulence models accurately predict integral quantities, mean velocity profiles and turbulence 
intensity profiles in the wake. Although some quantities are predicted slightly better by a particular model, 
the difference is not sufficient enough for one to assert that one of the models is the best. For subcritical 
Reynolds numbers, it is expected that all DES models should behave like a good LES model. Further 
numerical tests, specifically, for wall bounded internal flows where turbulent separation and reattachment 
locations are more difficult to predict, remain for future study. 
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